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Abstract. We discuss a family of models expressed by nonlinear differential equation systems describing 
closed market societies in the presence of taxation and redistribution. We focus in particular on three 
example models obtained in correspondence to different parameter choices. We analyse the influence of 
the various choices on the long time shape of the income distribution. Several simulations suggest that 
behavioral heterogeneity among the individuals plays a definite role in the formation of fat tails of the 
asymptotic stationary distributions. This is in agreement with results found with different approaches and 
techniques. We also show that an excellent fit for the computational outputs of our models is provided by 
the K-generalized distribution introduced by Kaniadakis in [14] . 

lenge is to derive from the knowledge of the microscopic 
interactions the emergence of collective patterns similar 
to those observed in the real world. The techniques more 
frequently employed take advantage of today's computer 
power and include Monte Carlo methods and agent based 
or numerical simulations. 

In a recent work [8j belonging to this research line the 
authors of this note investigated a family of models for 
the taxation and redistribution process, formulated within 
a general framework proposed in The models at is- 
sue are expressed by systems of nonlinear ordinary dif- 



1 Introduction 

Several papers have been published, especially during the 
last decade, which aim at explaining some universal fea- 
tures of income and wealth distributions with the tools 
of statistical mechanics and kinetic theory. We recall for 
example [J, g], [3J, @], [5], and refer to [5] or [7] for a 
detailed review and bibliography. The models described 
in most of these papers typically concern closed market 
societies in which wealth exchanges obeying some given 
rule and involving pairs of agents take place. The chal- 
Send offprint requests to: 
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ferential equations, which enjoy a great descriptive flexi- scribed by a (continuous) Boltzmann-like equation, which 
bility thanks to the presence of various parameters. The proves in particular that randomness in the microscopic 
equation systems are well beyond analytic solution. Never- trade and saving propensity play a significant role in the 
theless, many computational simulations give evidence of large wealth behavior of the steady distributions, 
the existence of asymptotic stationary distributions whose Taking advantage of the high flexibility which charac- 
higher income sections exhibit, at least for some choices terizes the models in [8], we give account below of several 
of the model parameters, a power law like decreasing be- numerical simulations run in correspondence to various 
havior. Together with other aspects discussed in [8] and choices of parameters. Modulo the different model formu- 
la, this fact encourages to pursue the study of these mod- lation, these simulations suggest, in substantial agreement 
els. Indeed, as well known after Pareto observation of more w ith the above mentioned references, that an heteroge- 
than a century ago [10 j, often in real economies the wealth neous saving propensity is a crucial ingredient for the fat- 
distribution of the richest people follows a power law pat- ne ss of the large wealth part of the stationary distribution. 



ternQ 



A second point in this paper concerns the attempt to 



In particular, the question arises of what is within the find the analytical expression of a multi-parameter func- 

framework and the models introduced in [8] the cause of tion with which to fit the equilibrium distributions ob- 

the detected decreasing behavior. This is the first question tained via numerical simulations. Various possibilities sug- 

addressed in this paper. A mechanism responsible for the gested in the literature are briefly recalled in Section [4] A 

appearence of a power law tail was recognized in [11] and comparison among them singles out the K-generalized dis- 

(U for a class of kinetic wealth exchange models to be the tribution introduced and discussed in [Hj , [TS] and [Hj as 

so called saving propensity, a parameter expressing the at- especially suitable to fit at a time the computational out- 

titude of agents to retain, when trading, a certain fraction puts relative to the poorest, average and richest income 



of their wealth. An analysis of the influence on the shape of values of the asymptotic distributions of our models 
the distribution of homogeneous and heterogeneous (i.e., 
equal or diversified) saving propensity among agents of a 
population can be found also in |12j . where various ki- 
netic wealth exchange models of both types are surveyed. 
Furthermore, in [13] a rigorous analytic investigation is 
carried out for a class of conservative economy models de- 



The plan of the paper is as follows. In Section[2]we re- 
call the general framework introduced in ]9 and discuss a 
family of models which enjoy minor changes with respect 
to the one in [8]. In Section|3]we particularize the models, 
by choosing three specific sets of values for their parame- 
ters and we focus on the effect produced by the different 



1 The words income and wealth are used in this paper to choices on the shape of the asymptotic income distribu- 
designate the same concept. tions. In the same section we also review and shortly com- 
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pare some other wealth exchange models discussed in the through taxation to provide welfare services for the popu- 
econophysics literature. In Section [4] we extend the atten- lation, we interpret the welfare provision as an income re- 
tion to the lower income parts of the asymptotic income distribution. Ignoring the passages to and from the govern- 
distributions found and we show, for each of the models ment, we adopt the following equivalent mechanism as the 
of the previous section, how the "numerical" distributions mover of the dynamics: in correspondence to any interac- 
fit with suitably tuned K-generalized distributions. Some tion between an /i-individual and a /c-individual, where the 
remarks on reversibility can be found in Section [5j In Sec- one who has to pay S to the other one is the /i-individual, 
tion [6] the results are summarized and discussed within a this pays to the /c-individual a quantity S (1 — r) and he 
critical perspective. pays as well a quantity St, which is divided among all j- 

individuals for j ^ nj^j Accordingly, the effect of taxation 
2 The general framework and a model family and redistribution is equivalent to the effect of a quantity 

of interactions between the /i-individual and each one of 



Consider a population of individuals divided into a finite 
number n of classes, each one characterized by its average 
income. Let ri, r 2 . . . r n denote the average incomes of the 
n classes, ordered so that r\ < r 2 < . . . < r n , and let 
Xi(t), where Xi : R — > [0, +oo) for i e F n = {1,2, —,n}, 
denote the fraction at time t of individuals belonging to 
the i-th class. In the following, the indices i,j,h,k, etc. 
always belong to r n , if no differently stated. 

Assume that pairwise interactions of economic nature, 
subjected to taxation, take place. And call S the fixed 
amount of money that people may exchange during their 
interactions. Vhk € [0, +oo) , 

Any time an individual of the h-th class has to pay a 
quantity S to an individual of the fc-th class, this one in 
turn has to pay some tax corresponding to a percentage of 



the j-individuals for j ^ n, which are "induced" by the 
effective h-k interaction. To fix notations, we may distin- 
guish between direct interactions (h-k) and indirect inter- 
actions (h-j for j ^ ri). 

Any direct or indirect economical interaction yields as 
a consequence a possible slight increase or slight decrease 
of the income of individuals. 

To translate into mathematical terms all that, we in- 
troduce 

- the (table of the) interaction rates 



expressing the number of effective encounters per unit 
time between individuals of the h-th class and individuals 
of the fc-th class; 



what he is receiving. This tax is quantified as S r, with the 2 The reason why individuals of the n-th class constitute an 

tax rate r = r k < 1 depending in general on the class of exception is a technical one: if an individual of the n-th class 

the earning individual. Since the quantity S r goes to the would receive some money, the possibility would arise for him 

government, which is supposed to use the money collected to advance to a higher class, which is impossible. 
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- the (tables of the) direct transition probability den- In order to design within the general framework a 

sities specific model (or a specific family of models), we need to 

C hk £ [0, +00) , further characterize the expressions of the direct transi- 

satisfying for any fixed h and k tion Probability densities C\ k and the indirect transition 



i=i 



variation densities T^ hk ^ (x) . A conceivable choice is given 
next. 



which express the probability density that an individual TTT . .. . . 

We represent the direct transition probability densities 

of the h-th class will belong to the i-th class after a direct ni 

Oll as 



interaction with an individual of the fc-th class; 

- the (tables of the) indirect transition variation den 
sities 

1 [hk] H ' H 

where the T? hk Ax) with x — (x\, ...,x n ) £ R™ are contin- 
uous functions, satisfying, for any fixed h, k and x € R n 



Chic — a hk + frftfc > 



where the term a hk expresses the probability density that 
an /i-individual will belong to the i-th class after an en- 
counter with a fc-individual, when such an encounter does 
not produce any change of class and the term b hk expresses 
the density variation in the i-th class of an /i-individual 

n 

interacting with a fc-individual. 

i=l 

„ ,/> ,t *j* 1 Accordingly, the only nonzero elements al h are 

I hese functions account lor the indirect interactions and nh 



express the variation density in the i-th class due to an 
interaction between an individual of the h-th class with 
an individual of the fc-th class. 

Chose for simplicity all the interactions rates rjhk to 
be equal to 1, which corresponds to assuming that all the 
encounters between two individuals occur with the same 
frequency, independently of the classes to which the two 
belong. 

Then, the evolution of the class populations Xi(t) is 
governed by the following differential equations, in which 
the contribution of both direct and indirect interactions is 
present: 

, n n n 



<% = 1 



To define the elements b l hk , we introduce the matrix 
P, whose elements ph,k express the probability that in an 
encounter between an /i-individual and a fc-individual, the 
one who pays is the ^.-individual. In view of the possibil- 
ity that to some extent the two individuals do not really 
interact, the pn,k are required to satisfy < pu tk < 1 and, 
of course, Ph,k + Pk.h < 1- Apart from that, there is a 
certain arbitrariness in the construction of the matrix P. 
The only additional requirement is to put each element 
in the first row, as well as each element but the very last 
one in the last column, equal to zero. The reason for that 



dt f— ' f— ' V / f— ' is the non existence of a class lower than the first nor a 

/l — 1 K — 1 fc — 1 
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class higher than the n-th: we cannot admit the possi- We express the indirect transition variation densities 
bility for 1-individuals [respectively, for n-individuals] to (x) as 

move back to a lower class [respectively, to advance, pass- , „ , , , „ 

T[ hk] (x) = U{ hk] (x) + V [hk] (x), 

ing in a higher class]. Observing that in the presence of 

where 

interactions between two individuals of the same class, the 

Ph,k S T k ( Xj_l 

r i+ i-ri 

into account the direct transition, and may only change 

represents the variation density corresponding to the ad- 

because of the taxation and redistribution contributions, 

vancement from a class to the subsequent one, due to the 

we assume that individuals of class 1 never pay. We also 



average wealth of the two remains the same if we only keep ^[hk] ( x ) = ~ ( ' ' ) (^) 



benefit of taxation and 



assume that individuals of class n never receive money 



V[hk]( x ) =Ph,kSr k 



An encounter between an ft-individual and a fc-indi- ' \ r h — fi Th — / 2j=i 

(4) 

vidual where h ^ k, with h > 2 and k < n — 1 and the 

_ with 5h,k denoting the Kronecker delta, accounts for the 

/i-individual paying, produces the elementfrl 

variation density corresponding to the retrocession from 

uh—l _ c ~ Tk h h — — o ~ Tfc 

°hk ~ Ph,k ° „ _ „ ' °hk ~ Ph.k 1 ) a class to the preceding one, due to the payment of some 

' h * h—1 ' h ' h—1 

bXX 1 =PhkS — — T — , bi h = -ph k S — — — . tax. In the r.h.s. of p| and pi, ft > 1 and the terms 
rk+i-r k " ' r k+1 -r k u ^ 

involving the index i — 1 [respectively, i + 1] are effectively 

Hence, the possibly nonzero elements b hk are of the form 

present only provided i — 1 > 1 [respectively, i + 1 < n] . 
, = fe 5 — — , Notice that for technical reasons, in the model under 

ii _ o 1 — T * _ c — Tfc consideration, the effective amount of money paid as tax 

b ik = ~Pk.i^ Pi,kd ) ' 

n+i -n n- n-i 

1 — Tj_i relative to an exchange of 5(1 — r k ) between two indi- 

°i-i,k = Pk,i-i S , (2) 

viduals and then redistributed among classes is given by 



where the expression for b\ +l k in (j2j) holds true for i + 1 + 5rfe ^ Xj )/^ =1 Xj ) instead of Sr k . 

k, and for i < n— 1 and k < n— 1; in the expression for b) k , . _ . 

A general theorem proved m [9j ensures that, with the 

which is only present whenever i ^ fc, the first addendum . 

present choice of parameters, m correspondence to any 

initial condition xo = (xqi, . . . , xo n ), for which xoi > and 

Y^i=o x oi = lj a unique solution x(t) = (xi(t), . . . ,x n (t)) 

and the expression for b}_i h holds true for for i — 1 ^ k, „ r-n . , . r „ \ 

of |1J) exists, which is defined for all t e [0, +oo), satisfies 



is effectively present only provided i < n— 1 and k > 2 and 
the second addendum only provided i > 2 and k < n — 1, 



and for i > 2 and fc > 2. 

3 Since the classes are characterized by an average income, 
the terms b\ h are equal to zero for any h and i. 



x(0) = xo and also 

n 

x t (t) > for i e r n and Xj(t) = 1 for all t > . (5) 
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n n 

Y Y iphk + T [hk] ( x )) x h x k -Xi, i£ r n , (6) 



With reference to this well-posedness result, we observe 
that in this context the only meaningful initial data, and 
solutions as well, are those with non negative components. 
Also, the constraint X)"=o x i0 ~ •"■ expresses nothing but 
a normalization. Hence, we conclude by ([5| that the solu- 
tions of interest are in fact distribution functions. Further- 
more, we emphasize that in view of the expressions of 
the U* hk ^ (x) and Vj^fe] ( x ) m Q and Q can be simplified. 

We may from now on consider, instead of ([lj, the sys- 
tem of differential equations 

dxj _ 
dt ~ 

h=l k=l 

where the terms Tf^ fc i (x) are linear in the variables Xj . 
The equations in ([6| have a polynomial right hand side, 
containing cubic terms as the highest degree ones. 

Due to the fact that the value of n and the parameters 
rk,Tk,Ph,k are still to be fixed, the equations |6]) actually 
describe a family of models rather than a single model. 

3 Different choices of the parameters 

We now exploit the flexibility which characterizes the equa- 
tions ^ and take into consideration three different possi- 
ble choices for the values of the parameters Ph,k- Our aim 
is to try and see what is the effect of such different choices 
on the shape of the large wealth part of the asymptotic 
income distribution. 

3.1 Three different cases 

In each of the three considered cases we take the elements 
of the matrix P lying on the main diagonal, those on the 



first row and those on the n-th column to be equal to zero: 

p jt j=0 for j £ {l,...,n}, 
pi tk = for k € {l,...,n}, 
Ph.n = for h € {1, ...,n} . 

Apart from that, we define the pn^k 

• in a first case, 

which we'll call the Case FSP, 
(FSP for fixed saving propensity), 
as 

Ph,k = (l/4)r fe /r„ 

with the exception of the terms 

Ph,l = (1/2) r h /r n for h € {2, n} , 
Pn.k = 1/2 for k £ {1, ...,n — 1} ; 

• in a second case, 

which we'll call the Case TID, 

(TID for trading individuals dependence), 

as 

PhM = (1/4) min{r h ,r k }/r n 

with the exception of the terms 

p h<1 = (1/2) n/r„ for h £ {2, n} , 
Pn,k = (1/2) r k /r n for k £ {1, n - 1} ; 

• in a third case, 

which we'll call the Case FAE, 
(FAE for fixed amount exchange) , 



as 



Ph.k = 1/4 
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with the exception of the terms The influence on such distribution of the different in- 
teraction rules relative to the Cases FSP. TID and FAE 

p M = l/2 for he{2,...,n}, 

is well illustrated by the Figures 1,2, referring to the case 

p n ,k = 1/2 for k e {1, ...,n - 1} ; 

n = 25, just a sample among several similar ones. In each 
Choosing the p h ^ k as in the Case FSP amounts to as- of these figures two panels relative to the Case FSP in 
sume, in view of that each individual (apart from the first row, two relative to the Case TID in the sec- 
those belonging to the poorest and the richest class, for ond row, and two relative to the Case FAE in the third 
which particular rules hold true) has the same saving pro- row are reported. The panels on the left show the his- 
pensity. Indeed, the frequency of interactions is uniform tograms corresponding to an initial condition; the panels 
and the effect of the mechanism described above is the on the right show the corresponding asymptotic station- 
same as if in any trade each individual would pay a small ary distribution. We point out here that the histograms 
quantity proportional to her income, the proportion rate are scaled differently from panel to panel. The initial con- 
being the same for everybody. Such a homogeneity in ditions are chosen randomly within distributions for which 
spending and hence in saving is replaced by heterogeneity lower income classes contain a high percentage of individ- 
both in the Cases TID and FAE. In particular, the Case uals, which seems to be a rather realistic assumption. In 
TID expresses the conjecture that, when two individuals each of the two figures they are taken to be the same 
trade, the exchanged amount of money is proportional to for the three simulations (each one referring to a differ- 
the wealth of the poorest of the two. ent model) precisely to emphasize the role played by the 
The average incomes r 6 and the taxations rates Tj are different P aramet ers in the formation of the shape of the 
taken here to be given by r 3 = Wj and asymptotic distribution. 



j _ i The Figures 1 and 2 suggest that a fat tail exists in 

Train) 

the Cases TID and FAE, and not in the Case FSP. As a 



n — 1 



for j e r„, where T min = 10/100 and T max = 45/100. m0 re, a log-log plot of the right section of the asymptotic 

» j. r i i r ■ stationary distributions is in agreement with a power law 

According to the result ot a large number ot simu- J ° 

, ,. / i I7T1 j 1 1 i r i ,i type decrease for the Cases TID and FAE, while this does 

lations (see also [9] and [8jJ, once the value of n and the Jt ^ 

n , . n not happen in the Case FSP. 
parameters Tj and Tj are chosen, m correspondence to any 

fixed value /x £ [r±, r n ] of the global wealth a stationary To gain some accuracy, we tried several simulations 

distribution exists, which is the asymptotic trend of all in correspondence to a different choice of the average in- 

solutions of ^ with initial conditions x a = (xqi, . . . , Xq tl ) comes: in particular, we took them to be nonlinear: Tj = 

satisfying x oi > 0, J2i=o x 0t = 1 an d Ya=i riXoi = A 4 - 10(1,67)- J for j G T 25 , so as to have r 2 4/ri 5 — 100. The 



8 
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Fig. 1. Initial (on the left) and long-time equilibrium (on the right) distributions for the models at issue: the first row 
refers to the Case FSP {fixed saving propensity), the second row to the Case TID [trading individuals dependence) , 
the third row to the Case FAE (fixed amount exchange). Notice that the histograms are scaled differently on different 
pictures. 




Fig. 2. Initial (on the left) and long-time equilibrium (on the right) distributions for the models at issue: the first row 
refers to the Case FSP (fixed saving propensity), the second row to the Case TID (trading individuals dependence), 
the third row to the Case FAE (fixed amount exchange). Notice that the histograms are scaled differently on different 
pictures. 
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expectation was confirmed: no power law behavior is seen computed explicitly from the simulations tends to coincide 
in the Case FSP, while a clear straight line appears in the for large k with the integral of the rescaled density, pro- 
log-log plots of the Cases TID and FAE. As a sample of vided one also takes into account certain slight differences 
the result of several simulations, three log-log plots rela- in the normalization (the discrete incomes are normalized 
tives to the Cases FSP, TID and FAE are shown in Figure as Y^i=i x i = 1> while the rescaled density f(r) is normal- 
3. ized as J °° f[r)dr = 1). 

In the simulations for the investigation described here 
3.2 Some remarks on finding power laws we employed, besides linear class incomes, incomes which 

increase exponentially as = 10 (1,67)*. The basis 1.67 

In general, simulations based on classes whose income val- 

is chosen in such a way as to obtain an income increase 

ues increase linearly do not allow a precise determination 

by a factor 100 when the class index i increases by 10 

of the tail of the equilibrium distribution. It is well known 

(the tail of the distribution is typically fitted to look for 

(see e.g. [17]) that in real statistical data a fraction of less 

a power law in the range from i = 15 to i = 25). The 

than 2-3 % of the population has an income distributed 

rescaling factor to be applied to the discrete distribution 

according to a power law. But a power law tail can be 

function Xi is in this case proportional to the reciprocal 

clearly distinguished from an exponential tail only over 

income ; this means that a power law tail maintains, 

an income range of two or three magnitude orders. This 

after rescaling, the same power law form, with an effective 

requires the introduction in our model of class incomes 

Pareto index equal to the index of the Xi augmented by 1. 

growing more steeply than linearly; in this way the indi- 
viduals of the 25-th class would not be, for instance, 25 

times richer than those of the first class, but at least 1000 3.3 A short survey of some wealth exchange models 

times richer. 

In our previous work [5] we introduced for this pur- To put our analysis into an historical prespective, we re- 
pose class incomes which increase quadratically; in that call here that several mathematical models of wealth ex- 
case the income of the 25-th class was 25 2 times larger change processes display a Pareto tail in their equilibrium 
than the income of the first class. Notice that when the solutions. In some cases the tail has only be observed nu- 
class income ri is not linear in the z'-th class, a rescaling of mcrically, in others it is a part of an analytical solution, 
the discrete distribution function Xi is necessary, in order In some models the Pareto index is a constant, in oth- 
to preserve the correct relation between the distribution ers it depends on some parameters. Finally, some models 
density and the cumulative income distribution [8], It is are conservative, the conserved quantity being the total 
possible to check that the cumulative distribution 53i=i x i wealth; in others the total wealth increases in time. 
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Fig. 3. The (suitably rescaled) log-log plots of the right section of asymptotic stationary distributions (in a case with 
average incomes rj = 10 (1, 67)-?) in the three Cases: FSP, TID, and FAE. A clear straight line appears in the log-log 
plots of the Cases TID and FAE, while this does not happen in the Case FSP. 

For instance, the model by Bouchaud and Mezard |18j . turn. One finds in this case a power law distribution for 

originally developed in physics for the description of poly- the rich of the form f{x) ~ a; -1 ' 7 . 

mers, admits (when the exchange rate J is the same for l n the classical conservative model by Chatterjee, Cha- 
all agents) an exact solution in the mean-field approxima- krabarti and Manna [11], agents contribute only a frac- 
tion. The distribution function is in that case an Amoroso tion of their wealth for trading, depending on their saving 
curve with a Pareto index a = 1 + J /a 2 which depends on propensities, which differ among agents. A Pareto tail is 



the exchange rate J and on the variance a of a Gaussian found, of the form f(x) ~ x~ 



2 



multiplicative process appearing in the model equations. 0ur model family som ehow resembles that one by Chat- 

The Gaussian process simulates investment dynamics and ter j ee; Chakrabarti and Manna [TT]. However, differently 

is related to the temporal change in the value of the stocks from itj it allows to recover different Pareto indices (com- 

which, on the average, increases in time. pare e g results in jg; or see the data in the caption of the 

The model by Scafetta, West and Picozzi [19] is an Figure 4). 
evolution of the model by Bouchaud and Mezard, which The question of taxation and redistribution was first 

takes into account several realistic details of the exchange studied in [T], where it was modeled through a Boltzmann 

process. Its computer-generated wealth distribution func- equation. A different approach was developed in [21], en- 

tions can be fitted by rational functions of the form f(x) = compassing two-step tradings, which consist of a wealth 

axt/ (1 + bx)'< +& , where the Pareto index is a = 5 - 1. exchange analogous to an inelastic binary "collision "and 

In the model by Slanina [20] , growing markets are mod- of a redistribution of the lost wealth (the taxes) among 

eled by bringing in during each trading an extra amount of the population. The effect of the subsidies by the gov- 

wealth, which is proportional to the wealth of both agents ernment on the equilibrium distribution is shown in both 

participating in the trading: this represents, realistically, papers pQ and |21) to cause a shifting of the individu- 

the fact that the rich can invest more, and with more re- als from the lower income classes toward middle income 
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classes. Accordingly, and differently with respect to what presence of a maximum, by a null value in zero (with ex- 
happens e.g. to Boltzmann-Gibbs distributions, the equi- ponential increase) and by a Pareto index a which is in- 
librium distribution is seen to exhibit a maximum. The dependent from the average income. 

Figure 5 of [T] and the Figure 2 of [3T] are qualitatively The Pareto-Gompertz distribution contains three in- 
similar to the figures of this paper with the histograms dependent parameters, is monotonically decreasing and 
representing the asymptotic equilibria. has a Pareto index which depends on the average income. 

Its analytical expression is more complex: it is given by a 

4 The K distribution as a reasonable fit combination of a double exponential of the form f(x) = 



When the equilibrium income distribution /(x) of a model 
is found through numerical methods, it is customary to 
characterize this function according to its general features 
(e.g., monotonicity, presence of a maximum, value at zero) 
and to compute some global statistical indices, as the vari- 
ance, the Gini index, etc.. Finally, one tries to fit f(x) with 
some known multi-parameter distribution. 

It is well known, for instance, that reversible kinetic 

= a/3x a - 1 exp K (-f3x a ) 

systems evolve towards a Gibbs function, independently ^/l + /3 2 k 2 x 2q ' 

from the details of the transition probabilities, while non- with the generalized exponential function 
reversible systems with fixed saving propensity evolve to- expJy) = (\A + n 2 y 2 + ny) 1/K 



followed by a power tail. It has been suc- 
cessfully employed e.g. in [3] for the description of real 
statistical data from Brazil. 

Rational functions of the form f(x) — ax 1 / ' (l + bx) 1+s 
are suggested by Scafetta, West and Picozzi in [T5] as ap- 
parently well fittings the computer-generated wealth dis- 
tributions of their trade-investment model. 
The re-generalized distribution, given by 



wards a Gamma function. For systems with a Pareto power 
law tail, four are, to our knowledge, the fit functions which 
have been proposed in the literature: the Amoroso distri- 
bution |22j . the Pareto-Gompertz distribution [3], some 
rational functions of the form /(x) = ax 1 / (\ + bx) 1+s [lQ , 
and the re-generalized distribution [TJ], [15] , [15] . 
The Amoroso distribution 



was introduced by Kaniadakis in [14] in a general context. 
In [15] and [16] it has been considered in connection with 
its application to the economic problem of the representa- 
tion of income distribution, with the parameters a, /3 > 
and re € [0,1). In that context, the variable x in Q is 
defined as x = z/(z), with z the absolute personal income 
and (z) its mean value. One of the interesting properties 

(a - l) Q e"^ 

J\ x ) — — p7 — \ i +a that this probability density function enjoys is its behavior 



r{a) x 1+c 



for large values of x; indeed, 



contains one parameter and is the exact limit of the model 

by Bouchaud and Mezard [15]. It is characterized by the /(») ~ - (2/3 k)~ 1/k x~ ( for x +oo . 



/v 



12 



M.L. Bertotti and G. Modanese: A Family of Taxation and Redistribution Models 




Fig. 4. In each of the two panel rows three fittings are illustrated of three asymptotic stationary distributions relative 
to the three Cases FSP, TID, and FAE of Section [3] with the ^-generalized distribution. The taxations rates Tj 
are different in the simulations relative to the first and to the second row. The horizontal variable is defined as 
x = z/(z), with z the absolute personal income and (z) its mean value. The vertical variable is rescaled so as to have 
a normalized distribution. The parameters a, j3 e k of the ^-generalized distribution ffify are as follows. In the first row: 
FSP: a = 2.9459, /3 = 0.7754, k = 0.4426; TID: a = 1.5113, /3 = 1.0574, k = 0.6156; FAE: a = 1.4601, (3 = 1.3008, 
K = 0.9952; in the second row: FSP: a = 2.8868, (3 = 0.7776, k = 0.4342; TID: a = 1.5118, j3 = 1.0560, K = 0.6138; 
FAE: a = 1.3784, (3 = 1.2171, k = 0.8664. 



The parameters a and /3 respectively determine the cur- 
vature (shape) and the scale of the probability distribu- 
tion. The parameter k measures the fatness of the up- 
per tail: the larger its magnitude is, the fatter the tail 
is. The strength of the ^-generalized distribution stems 
from its ability to reproduce in an excellent way empiri- 
cal data of several countries (see [T5] for data concerning 
Germany, Italy, United Kingdom and |16j for Australia, 
United States). 

Based on several attempts to fit the asymptotic sta- 
tionary distributions resulting from our simulations with 
one or the other of the mentioned analytical expressions, 



we are in the position of claiming that the unique among 
them which seems to be suitable for the purpose is the 
K-generalized distribution. In fact, it displays quite good 
performances. 



The Figure 4 illustrates in each of the two panel rows a 
sample of three fittings out of many, relative to the three 
Cases FSP, TID and FAE introduced in Section [3] The 
average incomes are Tj = lOj. The taxations rates Tj are 
taken differently in the simulations relative to the first and 
to the second row. Specifically, in the simulations whose 
asymptotic trend is shown in the first row, three taxation 
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rates are postulated, given by: 



Tj = 5/100 for j = l,...,5, \.. < ....^.. S . < .... 

Tj = 25/100 for j = 6, ...,20, 



73=45/100 for j = 21, ...,25, 2 2 

while in the simulations relative to the second row the 



taxation rates are as in In each of the plots the hori- \ 
zontal variable is given by x = z/(z) , with z the absolute 
personal income and (z) its mean value and the vertical 
variable representing the fraction of individual with a cer- 
tain income is rescaled so as to have a distribution. 



/ 1 
/ 1 
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5 Some remarks on reversibility 



Fig. 5. Example of non-reversible process with taxes. The 
first diagram shows an example of an interaction between 
a "rich" agent (belonging to the 20-th income class) and 



According to general results of kinetic theory [B] any sys- a "poor" agent (6-th class). In this case the poor agent 

tern with reversible microscopic interactions approaches, P a y s an amount S and acquires a definite probability to 

at equilibrium, a Gibbs distribution function (a pure expo- P ass to the lower class - The rich a S ent S ains an a mount S 

nential). The reversibility condition is satisfied when the minus the taxes and acquires a probability to pass to the 

probability of each inverse elementary interaction process hi g her class - Thc taxes are distributed among all the other 

in the system is equal to the probability of the correspond- a S ents ' S ivin g each of them a ( sma11 ) probability of transi- 

ing direct process ^ on to their upper classes. The second diagram represents 

For our model family, a typical example of elementary the inverse process, where the rich agent pays the poor 

interaction with money exchange and taxation is shown a S ent and thus acquires a probability to pass to the lower 

in Fig. 5 (see detailed discussion below). This interaction class - Suppose that the two alternatives "6 pays 20" and 

is reversible in only one case among those described in " 21 P avs 5 "have tne same probability. If the rich and the 

this paper, and more precisely only provided thc following P 00r are subjected to different taxation rates, then the two 

three conditions are satisfied (see Table 1): diagrams have different probabilities and the whole pro- 
cess is clearly non-reversible. But even if the taxation rate 



1. the parameters fry are as in the Case FAEoi the model 
family, in which the exchanged amount is fixed and the 
saving propensity is variable; 

2. the class income values T\- are linear in the class in- 
dex j; 



is the same, the whole process is non-reversible, because 
there exists no way for all the other agents to give the tax 
money back. In other words, the exchange of the amount 
S can be reversed, but not the "small re-distributions" of 
the tax on S. 
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3. the taxation rates are set equal to zero. agents, giving each of them a (small) probability of transi- 

tion to their upper classes. The second diagram represents 



Let us analyse these three conditions in detail. 

As for condition 1, it is well known |6] and easily re- 
alised that wealth exchange processes with fixed saving 
propensity are non-reversible. Such processes are present 
in the Case FSPof our model, where phk = (1/4) Thjr n . In 
Case TID, where phk = (1/4) minjr^/rfel/rn, the interac- 
tion is also non-reversible, because the quantity min{r/j/rfc} 
is in general different in the direct process and in the in- 
verse process. 

In connection with condition 2, we notice that the in- 
come differences |rj — r^_i| in the denominators of the pa- 
rameters 6^ in ([2]) are independent from the index i only 
if Ti is linear in i. Otherwise, the income differences de- 
pend on income, being typically larger for the rich classes. 
It follows that also in that case the probabilities of the 
direct and inverse wealth exchange processes are not the 
same. 



the inverse process, where the rich agent pays the poor 
agent and thus acquires a probability to pass to the lower 
class. Plainly, if the rich and the poor are subjected to dif- 
ferent taxation rates, then the two diagrams have different 
probabilities and the whole process is non-reversible. But 
even if the taxation rate is the same, the whole process 
is non-reversible, because there exists no way for all the 
other agents to give the tax money back. In other words, 
the exchange of the amount S can be reversed, but not 
the "small redistributions" of the tax on S. 

Hence, calculating the distribution function of our mo- 
del in Case FAE, with linear income classes and no taxes, 
provides an important cross-check of the model and of the 
numerical solutions: we expect in this case to obtain the 
pure Gibbs distribution predicted by kinetic theory for 
reversible processes. And this is, actually, what is found: 
one obtains an exponential distribution, except for the 
first and last classes, whose wealth exchanges are non- 
reversible due to boundary effects. 



Finally, about condition 3, even assuming that the con- 
ditions conditions 1 and 2 hold true, we have irreversibility 
due to the taxation process. This can be grasped by look- 
ing at the Figure 5, which is supposed to refer to Case 
FAE of our model family, with linear class incomes. The In this connection we also recall that the initial condi- 
first diagram shows an example of an interaction between tions play in general a role in the equilibrium approach to 
a "rich" agent (belonging to the 20-th income class) and a Gibbs distribution: such an equilibrium distribution can 
a "poor" agent (6-th class) . In this case the poor agent be reached in a conservative system only if the initial con- 
pays an amount S and acquires a definite probability to ditions do not display a population inversion correspond- 
pass to the lower class. The rich agent gains an amount S ing to negative temperature. In our case this is guaranteed 
minus the taxes and acquires a probability to pass to the by the choice of an initial population concentrated in the 
higher class. The taxes are distributed among all the other classes with low income (compare Section [3]). 
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The emergence, in the reversible case without taxes, of 
a purely exponential distribution is at the same time nat- 
ural and puzzling. One can naturally expect, in fact, that 
the removal of taxes causes an increase in the population 
of the poorest at the detriment of the middle classes, as 
it happens in the exponential distribution. On the other 
hand, the exponential tail is thinner than the power law 
tail of the (non-reversible) distribution with taxes, and the 
reason for this is less obvious. Things go as if the complete 
reversibility of the wealth exchange would make it harder 
for the super-rich to preserve their wealth. 

We close this section recalling that we pointed out 
in Section [3] a similarity between the Case FAE (with a 
fixed exchanged amount of money) of our model family 
and some models without taxes and with variable saving 
propensity ([H], see also [12] and [13]). We observed that 
our fixed exchanged amount of money mimics a saving 
propensity which varies among agents belonging to differ- 
ent classes and that in both cases a Pareto tail is found. 
Since, however, the removal of the taxes from our model 
causes the Pareto tail to disappear, one should conclude 
that the dynamical mechanism behind the formation of 
the tail is different in the two cases. 



Table 1. Summary of results. The three basic model cases 
considered in this paper (first column) and their further 
subdivision/specialization according to taxation and lin- 
ear/nonlinear class incomes. The last two columns record 
the appearence of reversibility and power law tails. Notice 
that linear class incomes in general do not allow to distin- 
guish clearly between a power law tail and an exponential 
tail. 



Model 


Income 


Reversibility 


Fat tail 




linear 


no 


? 


FSP with taxes 


non-lin. 


no 


no 




linear 


no 


? 


FSP, no taxes 


non-lin. 


no 


no 




linear 


110 


? 


TID with taxes 


non-lin. 


110 


yes 




linear 


no 


? 


TID, no taxes 


non-lin. 


no 


yes 




linear 


110 


? 


FAE with taxes 


non-lin. 


no 


yes 




linear 


yes 


no 


FAE, no taxes 


non-lin. 


no 


yes 



6 Concluding remarks 



exploited for the construction of various examples. A com- 
parison between different long time behaviors exhibited 
The aim of this paper is to provide further insight into the by those examples contributes to further understanding 



formation process of stationary income profiles. A family of the mechanisms responsible for the emergence of one or 

of models introduced in [3] and [5] has been investigated, the other output. Attention has also been given to plac- 

In particular, the high flexibility of these models (related ing our results within a frame of related results obtained 

to the wide range of possible parameter choices) has been through different methods. 
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Our approach for the description of the process of statements concerning the power law decreasing behavior 
wealth exchange in the presence of taxation comprises a are not to be intended in a strictly analytical way. We 
subdivision of the population into income classes, and a rather observe that for a more credible depiction of mod- 
corresponding discretization of the distribution function, ern market economies, more relevant corrections should 
This approach has proven to be powerful in its application, be incorporated in the model, the first one being proba- 
It translates well-tested physical methods of kinetic the- bly the introduction of a term accounting for a possible 
ory into a manageable mathematical framework (a system generation of wealth due to production efficiency increase, 
of ordinary differential equations) which can be quickly In fact, the conservation of the total wealth is not a real- 
solved numerically with standard software. The micro- istic assumption. 

scopic interaction parameters as well as the taxation rates We are very far from thinking and claiming that the 

have a straightforward interpretation and can be easily models investigated here can be employed to explain in a 

varied. This allows to reproduce the qualitative results of satisfactory way real world situations. However, they may 

several different models. provide a starting point for further progress in this direc- 
tion. The reason for their interest, we think, is given by 



In particular, in correspondence to some parameter 
choices, stationary distributions have been found, whose 
higher income sections exhibit a power law decreasing be- 
havior. It could be argued in this connection that properly 
talking of power law tails requires treating the income 
or wealth variable as a continuous one taking values in 
[0, +00). This is for example what is done in papers like 
[5 : and 1 - jl , where a rigorous analytical approach is devel- 
oped, which includes the study of the evolution of the mo- 
ments of the solution. Certainly, admitting a finite number 
of income classes entails some unnatural constraint and 
causes some boundary effect as well. On the other hand, 



their flexibility and ability of reproducing the emergence 
of qualitative global features based only on the description 
of the micro-level interactions. And this is one of the main 
challenges in a complex systems approach. 
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